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Abstract 



We consider the production, at future lepton colliders, of final fermion, 
sfermion, scalar pairs in SUSY models. For third family fermions and 
sfermions and for charged Higgses, the leading Yukawa effect at one loop for 
5^ I large cm. energies comes from a linear logarithm of Sudakov type, that only 

depends, in the MSSM, on one SUSY mass scale and on tan /?. Assuming 
a relatively light SUSY scenario, we illustrate a possible determination of 
tan P at cm. energies of about 1 TeV, working systematically at subleading 
logarithmic accuracy, at the one-loop level. 
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I. INTRODUCTION 



The possibility that virtual electroweak effects of Supersymmetric models at future 
lepton colliders [1,2] can be described, for sufficiently large cm. energies, by a logarithmic 
expansion of Sudakov type, has been examined in recent papers [3-5] in the simplest case 
of the MSSM. The first considered process has been that of final fermion pair production. 
This has been only treated at one-loop, computing for both massless [3] and massive [4] 
final states the SUSY Sudakov terms. One main conclusion is that these terms do exist 
and are all of "subleading" SL (hnear) kind. At one TeV, working in a relatively "hght" 
SUSY scenario, where the heaviest SUSY mass if of a few (two, three) hundred GeV, 
the numerical effects of SUSY Sudakov terms on the various observables are relatively 
small (a few percent at most), while at three TeV they become definitely large. The next 
step has been the analysis of final scalar (sfermion or Higgs) pair production for both 
massless and massive final states. This has been performed both at the one-loop level 
and at resummed subleading order accuracy [5]. By comparison of the two calculations 
it can be concluded that, in the assumed "light" SUSY scenario the two approximations 
are practically indistinguishable at one TeV, where their effect is, as in the fermion case, 
relatively small, but become drastically different (and both large, well beyond a relative 
ten percent) in the higher (two, three TeV) energy region. Within this approach one 
would thus conclude that at the LC extreme energies the MSSM can be safely treated, to 
subleading logarithmic order accuracy, at the one loop level for what concerns fermion, 
sfermion, scalar Higgs production. 

This conclusion could be of immediate practical consequence. In fact, it has been 
remarked in Refs. [4,5] that, for production of fermions and squarks of the third family and 
of charged Higgs bosons, the coefficients of the SL electroweak SUSY logarithms of Yukawa 
origin only depend on a common SUSY scale Ms, by definition the heaviest SUSY mass 
involved in the electroweak component of the process, and on the ratio of the two scalar 
vevs tan/9 = V2/V1. To the extent that a subleading order approximation can be considered 
as a reliable description of the variation with energy of the observables of the process, 
i.e. that the missing terms of the expansion can be adequately described by a constant 
component, this has allowed to propose a determination of tan/3 based on a number of 
measurements of the observables at different energies (roughly, on measurements of their 
slopes), whose main features have been already illustrated in Refs. [4,5] in a qualitative 
way for fermion production and scalar production separately. 

The aim of this note is that of proposing a more quantitative determination of tan/3 
from a combined analysis of the slopes in energy of fermion, sfermion, charged Higgs 
production. This will be done working at the one-loop level, in an energy region around 
(below) 1 TeV. With this purpose, and to try to be reasonably self-consistent, we briefly 
recall the structure of the various Yukawa contributions in the production of pairs of 
fermions (t,6), sfermions {tL,R, ^L,i?) and charged Higgs in the MSSM. The complete 
expressions of the asymptotic contributions can be found in Refs. [4,5], and we do not 
reproduce them here. Starting from them it is relatively straightforward to derive the 
quantities that are relevant for this note, which are given in the following list. 
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II. COMPLETE LIST OF YUKAWA EFFECTS IN CROSS SECTIONS AND 

ASYMMETRIES 



We parametrize the Yukawa effects in the physical observables that we are going to 
analyze and summarize them by giving a complete list. 

Let us denote by On, the various cross sections for production of sfermions {ti.R, ^Ljr), 
charged Higgs bosons and third generation fermions (t and b) . For top and bottom pro- 
duction we also include three basic asymmetry observables (unpolarized forward-backward 
asymmetry ApB, Left-Right asymmetry for longitudinally polarized beams A^r and 
its forward-backward asymmetry Ap^i). In the case of top production the average helicity 
Ht as well as its forward-backward and Left-Right asymmetries HFB,t and Hlr^i should 
be measurable by studying the leading top decay mode t Wh. The definition of these 
observables, in particular of the asymmetries, is conventional and can be found in full 
details in Appendix B of [6]. 

For cross sections O = cr, we define the relative one loop SUSY effect as the ratio 

O (n'^\ - COBorn+SM/ 2\ 
''n{q ) - 0Born+SM(^g2^ ' V^""^'' 

where "SM" denotes all the one loop terms that do not involve virtual SUSY partners 
(sfermions, gauginos and extra Higgs particles). For asymmetries, we consider instead the 
absolute SUSY effect defined as the difference 

en{q')^Oniq')-0^"'^'^^'^{q'). (2.2) 

At one loop, in the asymptotic regime, the shifts e„ can be parametrized as 

e(g^) = £ F(tan/3) In ^ + G + , (2.3) 

where, as we wrote, F is a simple function of tan f3 only. Its explicit expression must be 
determined by performing a Sudakov (logarithmic) expansion of the one loop calculation. 
The detailed analysis can be found in [3-5] and here we collect the various results for 
convenience of the reader. 



Sfermion cross sections 



F(^tJ = H cot' (3 + ml tan' /?j , (2.4) 

^K) = -yJ^ ml cot' (3, (2.5) 



F{a~J = F{ajJ, (2.6) 

,2 j.„„2 



n^J--lJ^rn-tta^'/3 (2.7) 
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Chcirged Higgs cross section 

F{aH±) = -jTT^ (m^t cot' + ^1 tan' /?) . (2.8) 



Fermion cross section 

(-3mi'(3 - 4s^ + 56s^) cot' (3 - ml{9 - 12s^ + 8s^) tan^ 

^^"^^ = .'^M^ 9 -24.^^ + 40.^^ ^ ^2-'°^ 
(-mi'(9 - 24s^ + 20s^) cot' /3 - 3m^(3 - 8s^ + 20s^) tan' /?) 



Fermion cross section asymmetries 



F{AFB,t) 
FiALR,t) 

F{Apoi,t) 

F{ALR,b) 

F{Apoi^h) 



m' cot' P — ml tan' P) 

1 384s^(3 — 4s^ + S^) , 2 ,2/3 2 J. 2 Q\ 

1 120.^^(9 - 12.;, + 8.^) _ 
(9 - 124, + 88st,)2 ^ 1 A' 6 

— 18s^/(3 — 8s^) ,2 ,2 ,3 2 . 2 ,3\ 
(9 -24.^, + 40.^^)2 ("^^ cot /? - tan /?) 

-96g^(3-8s^ + 5^^) 2 ^^^2 ^ _ ^2 ^^^2 ^) 
(9-24si^, + 40s^)2 ^ * ^ 
1 -305'^(9-24g^ + 205^) ^ 2 ,2^ 2, 



1 



Ml 
1 



(2.11) 
(2.12) 
(2.13) 
(2.14) 
(2.15) 
(2.16) 



Fermion helicity and its asymmetries 

It has been shown in [4] that the logarithmic parts of these observables are related to 
those of the cross section asymmetries as follows: 



F{Hf) = 




(2.17) 


F{Hfbj) = 


3 ^ 


(2.18) 


F{Hlrj) = 


4 , 


(2.19) 
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III. LIMITS AND CONFIDENCE REGIONS FOR TAN /? 



The constant G in Eq. (2.3) is a sub-subleading correction that does not increase with 
and depends on all mass ratios of virtual particles. The omitted terms in Eq. (2.3) 
vanish in the high energy limit [7]. 

To eliminate G we assume that a set of N independent measurements is available 
at cm. energies \f(^, • • • ? ^/^^ and take the difference of each measurement with 
respect to the one at lowest energy. For each observable, the resulting quantities 

Si^e{q^)-e{ql), (3.1) 

do not contain the constant term G and take the simple form 

di = F(tan/5*) In^, (3.2) 
Qi 

where tan /3* is the true unknown value that describes the experimental measurements. 

We now turn to a description of a possible strategy for the determination of tan f3. It 
results from a non linear analysis of data that must deal with extreme situations where 
tan (3 is determined with a still reasonable but rather large relative error. 

Let us label the various observables by the index n — 1, . . . , Nq and denote by cr„(g^) 
the experimental error on e„(g'^). For each set of explicit measurements {5n{Qi)}-i the best 
estimate for tan (3 is the value that minimizes the sum 

N No [i7;(tan/5)ln4i 
X^(tan/3) = E E , . (3-3) 

i=l n=l ^^n,i 

where = 5n{q^) and (T„^j = (Jn{q^). The factor 4 in the above formula follows from the 
fact that we assume a conservative error 2cr„ j on the difference In other words, we 
describe the experimentally measured quantity 5n,i in terms of a normal Gaussian random 
variable distributed around the theoretical value computed at tan /3* 

5n,i = K(tan/3*) In % + 2a„,, en,^, (3.4) 

with probability density for the independent fluctuations {^n,i} given by 

mn,i\) = Y{^e-'^^h. (3.5) 



In the following we shall simplify the analysis by taking a constant On^i = o with typical 
values around 1%. For each set of measurements we determine the optimal tan/3 that 
minimizes x^. It is a function of the actual measurements {Cn,i} and the width of its prob- 
ability distribution P(tan/3) determines the limits that can be assigned to the estimate 
of the unknown tan 
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The distribution P{tanf3) cannot be computed analytically because of the highly non 
linear dependence of the MSSM effects on tan/?. However, it can be easily obtained by 
Monte Carlo sampling. With this aim, we generate a large set of independent realizations 
of the measurements {Cn,i} and compute for each of them tan (3. The histogram of the 
obtained values is a numerical estimate of the true P(tan/3). 

In previous papers [4,5], we discussed a simplified approximate procedure and we 
determined the la boundary on tan j3 by linearizing the dependence on tan /3 around the 
minimum of x^. The bound that we derived is thus 

/No \"'^'/ o2 \-l/2 

5tan/3 = 2al^^F^(tan/3*)M (EI^'^J ■ (3-6) 

This result can be trusted if the experimental accuracy cr is small enough to determine 
a region around the minimum of where deviations from linearity can be neglected. It 
gives anyhow a rough idea of the easy regions where a determination of tan /3 from virtual 
one loop MSSM effects is not difficult. 

In a more realistic analysis, however, this approximation can be misleading and pos- 
sibly too much optimistic, especially for values of tan/? around 15 where the linearized 
analysis predicts typical relative errors around 50%. For this reasons, we pursue in this 
paper the complete Monte Carlo analysis of the allowed range of tan /3. 

IV. NUMERICAL ANALYSIS 

We begin by considering the full set of 16 observables consisting in: 

1. cross sections for sfermion production in the case of final ti, tR, bi, bL', 

2. cross section for charged Higgs production; 

3. cross sections and 3 asymmetries (forward-backward, longitudinal and polarized) 
for top, bottom production; 

4. 3 top helicity distributions (again forward-backward, longitudinal and polarized). 

We assume a set of = 10 measurements at energies between 600 GcV and 1 TeV 
with an aimed experimental precision equal to 1% for all observables at all energies. 
Within this ideal framework we have determined the probability distribution P(tan/3) 
for 5 < tan/5* < 40. In Fig. (1) we show the associated histograms in the four cases 
tan f5* = 10, 15, 20, 25. Even in the easiest case, tan f3* = 25, it is not possible to determine 
tan/? in a reasonable way. There is always a rather pronounced peak at small tan/3 in 
the histogram and the distribution is rather broad without a second peak recognizable 
around the exact value. 

To analyze in a quantitative way these results, we compute from each histogram the 
standard deviation of the estimated tan/3. If the distribution can be characterized in 
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terms of a single dominant peak, then this is rough measure of its width. Of course, when 
the distribution is wide or when it is the sum of two large separated peaks, the standard 
deviation is a pessimistic estimate of the uncertainty on the parameter determination 
that can be improved by adding, for instance, some information excluding the regions 
corresponding to large (or small) values. 

The determination of tan/? is almost completely driven by the observables related 
to sfermions and charged Higgses production. In Fig. (2), we show the results obtained 
without observables related to top and bottom production. In Fig. (3), we show the results 
obtained with top, bottom and charged Higgs observables. The single measurement of 
charged Higgs cross section is not enough to determine tan (5 with this level of precision 
in the measures. Finally, in Fig. (4), we show the results obtained with sfermions and 
fermions observables. The single measurement of charged Higgs cross section is not enough 
to determine tan j3 with this level of precision in the measurements. The plot of the relative 
error as a function of tan /3 for the various sets of observables (including the full case) is 
collected in Fig. (5). 

With these necessary remarks in mind, we can analyze the standard deviation of the 
parameter histograms and the result is shown in Fig. (5) (solid line). We see that for 
tan /5 < 20 a determination with a relative error smaller than 50% is not possible. 

If we consider still 10 measurements ranging from 600 GeV up to 1 TeV, but with a 
precision of 0.5%, then the scenario is quite better. In Fig. (6), we see that for tan /3* > 20, 
a well defined peak is visible in the rightmost part of the Figures roughly centered on the 
exact value. Fig. (10) (solid line) allows to conclude that the relative error is smaller 
than 50% for tan/?* > 13 and smaller than 25% for tan/?* > 25. Again, the role of the 
observables related to sparticle production is fundamental. In Fig. (7) we show what can 
be obtained without the information coming from top and bottom production. As one 
can see, there is small difference with respect to the previous two figures. Fig. (8) shows 
the results obtained with top, bottom and charged Higgs observables. Fig. (9) shows the 
results obtained with sfermions and fermions observables. Again, Fig. (10) collects the 
error as a function of tan /? for the various considered cases. 

As a final comment, we observe that a general feature of the histograms is the presence 
of a fake peak at small tan /? as well at tan /? ~ 6. The reason for this can be understood 
by analyzing what happens by exploiting in the analysis just the (dominant) charged 
Higgs cross section as discussed in Appendix A. 

V. CONCLUSIONS 

In a "standard" SUSY model, all the gauge couplings are fixed and coincide with the 
corresponding SM ones. For the couplings of the Yukawa sector, much more freedom is 
allowed. In the MSSM, one such coupling is the ratio of the scalar vevs tan (3. We have 
shown in this Note that, in a light SUSY scenario, a determination of tan/? based on 
measurements of the slope with energy of the combined set of observables of the three 
processes of fermion, sfermion, scalar charged Higgs production can lead to a determina- 
tion of this parameter with a relative error of 20-30 % in a range of high values that would 
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otherwise require difficult final state analyses of Higgs decays, see the proposals in [8]. 

The main point of our approach is, in our opinion, the fact that in this determination 
tan /3 is the only SUSY parameter to be measured: all the other parameters give vanish- 
ing contributions in the high energy limit. Isolating the various SUSY parameters to be 
studied is in fact, in our opinion, a basic feature of any realistic "determination strategy". 

In addition to the previous conclusion, we would like to add an extra final comment. If 
a SUSY model were different from the considered MSSM, in particular if it had a different 
Higgs structure (for example more Higgs doublets) , the Yukawa couplings would be, quite 
generally, different. But the features of the Sudakov structure would remain essentially 
unchanged. This would lead to the possibility of deriving, with minor changes in our 
approach, the components of the SUSY Yukawa sector that dominate the high energy be- 
haviour in this model. In analogy with what was done at LEPl for the "prediction" , from 
an analysis of one-loop effects, of the value of the top mass, that remains in our opinion 
one of the biggest LEPl achievements, a combined set of high precision measurements 
at future linear colliders physics could therefore produce a genuine determination of this 
fundamental SUSY parameter. 



APPENDIX A: ANALYSIS WITH THE CROSS SECTION ALONE 



It is interesting to analyze the role of charged Higgs production in tan /3 determination 
when no additional observal)los arc exi)l()ited. In fact, at the level of precision we are 
working, some problems arise. To see why, let us denote tan j3 = T and write explicitly: 

N [(F(T) - F(T*)) In ^ - 2aCi]'^ 



The function F{T) is given by 



PiT) = - ^ . {ml cot' (3 + ml tan' (5) (A2) 



If we denote Lj = log %i , the derivative of x' vanish when 

9i 



F'{T) = ^ T = ^mtjmi, ~ 6.2, (A3) 
and also at the solutions (if there are any) of 



F(T) - F(T*) = 2a^^^, (A4) 



that we can write in a simpler way in terms of a new normalized gaussian random variable 
I 



8 



A(r,r*) = ae (A5) 

where 

A(T, T*) = F(T) - F(T*), a = -j^^, (A6) 

To discuss the solutions of Eq. (A5), we must consider the main features of A(T, T*) for 
a given (unknown) T*. It tends to — oo for T ^ or oo and vanishes at 

Ti = T\ T2 = ^i^ (A7) 



It is a convex function and attains its maximum value aX T = rrit / rri}) where 

Each would-be measurement corresponds to a value of ^ and to an associated random 
T(^) that is found by minimization of x^- It is easy to see that two possibilities can arise: 

a) < l^maxiT*): in this case, X^{T) has a double well shape with a local maximum 



at T = ymt/mi, and two local minima around two points that are located around 
Ti and T2 and that tend toward them as a — 0. 

b) > Amax{T*)- in this case, X^(T) is concave and has a global minimum at T = 
^Jmt/mb. 

If the would-be measurements are randomly generated, cases (a) and (b) will occur with a 
relative frequency depending on a. For small a/Amax the majority of cases will be (a) and 
we shall be able to identify two well defined peaks in the histogram of the reconstructed 
T. The first will be false and around T2, the second will be true and around Ti = T*. Of 
course, if several measurements with independent dependencies on T are combined, then 
it is possible to suppress the false peak. 

If, on the other hand, a j A^ax is not small, then we shall fall in case (b) with very 
high probability and the reconstruction process will simply accumulate artificially at T = 



mt/rrih just because Eq. (A5) has no solutions. 
To give numerical values, with 10 measurements at 0.5% between 600 GeV and 1 TeV, 
we find that the condition a < Amax forbids the analysis of the region 3 < tan /3 < 13 
and in practice some other observable must be added (in the previous analysis we chose 
production of top or bottom) . 
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FIGURES 




FIG. 1. Histogram for the estimated tan^ minimizing x^- The experimental data included 
in the fit are all the observables discussed in the main text. We assume a precision cr = 1% on 
all data. The energy range is 0.6 < \/s < 1.0 TcV in this and the following figures. The four 
boxes show what happens at the particular values 10, 15, 20 and 25 of the parameter tan/3* 
that, we recall, is the true value. 
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FIG. 5. This figure shows the relative error in the estimate of tan/3 as a function of the true 
unknown value tan /?*. The four curves correspond to the different set of observables considered 
in the previous figures, as described in the legend. 
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